Quantum measurement theory for particle oscillations 
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A fundamental principle of quantum theory, clearly manifested in the two-slit experiment, is 
that for any alternatives that cannot be distinguished by measurement physical predictions are 
obtained by summation of their amplitudes. In particle oscillation experiments, a particle's time 
of detection is not directly measured, consequently, the detection probability should involve the 
summation over amplitudes corresponding to different detection times. However, in contrast to the 
principle above, standard treatments involve summation over probabilities rather than amplitudes; 
this implicitly assumes the existence of a decohering mechanism. In this work, we construct the 
detection probabilities for particle oscillations by summation over amplitudes, corresponding to 
different detection times. The resulting wavelength of particle oscillations differs from the standard 
expression by a factor of two. Moreover, we predict a dependence of the oscillation wavelength on 
the threshold of the decay process used for detection. 

PACS numbers: 03.65.Ta, 14.60.Pq, 03.65.Xp 
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The two-slit and related experiments demonstrate a 
striking feature of quantum mechanics: for any two al- 
ternatives that cannot be distinguished by a measure- 
ment scheme, their joint probability is obtained by sum- 
ming their contributions at the level of probability ampli- 
tudes; therefore interference terms appear. If, however, 
a macroscopic distinction of these alternatives is possi- 
ble, the interference terms are suppressed and the joint 
probability is the sum of the individual probabilities. 

Let us examine the importance of the above elementary 
point for the correct quantum treatment of particle oscil- 
lations [l[ . The detection time of the oscillating particles 
is not directly measured. Instead, we measure the num- 
ber of detection events, at distance L from the source. 
The question then arises: how, starting from the proba- 
bility amplitude A a {t, x) — for the detection of a flavor a, 
at a point x, at sharply defined moment of time t — , we 
can obtain the probability p a {L) that a particle of flavor 
a will be detected at distance L from the source. 

A common approach is to substitute in A a (t,x) the 
time t with the classical time of arrival at L, namely L/c 
for neutrinos, and to write the corresponding probability 
as p a (L) ~ \A a (L, L)\ 2 . The usual setting for this ap- 
proach is the textbook plane- wave treatment of neutrino 
oscillations, where 



Aa(t,x) = Up t U ai e 



a ipiX — iEit 



(1) 



Here, % labels the mass eigenstates of neutrinos, J7,- Q is the 
mixing matrix, (3 is the initial flavor, pi are the momenta 
in the different mass-cigenspaces and the corresponding 
energies E% = \/m 2 + p\ . The substitution t = L in Eq. 
(JTJ and the evaluation of the detection probability yields 
the standard expression for the oscillation wavenumbcr 
k i3 = (mf - mf )/(2p), for p t =p + 0(rrii). 



This method involves an unsatisfactory mixing of clas- 
sical and quantum concepts and this problem persists 
if we consider realistic wave packets rather than plane 
waves. More importantly, the method is ambiguous 
in the following way: different mass eigenstates corre- 
spond to different velocities Uj = Pi/Ei, if we evaluate 
each mass-component at the corresponding arrival times 
tj = L/vt, then we obtain an oscillation wave- number 
kij = (mf — Wj)/p for Pi = p + 0(rrii), which is twice the 
standard result 11 1 . 

Furthermore, the standard method assumes that the 
detection time should be computed in terms of an 'aver- 
age' velocity for all components of the wave packet — see, 
for example, Refs. [J, [5|. This assumption involves the 
substitution of a coherent superposition by an altogether 
different state that cannot be done in an invariant way 
in quantum mechanics. Any choice of an 'average veloc- 
ity' for a wave packet is arbitrary and unjustified at a 
fundamental level. 

An improved derivation of the standard expression for 
the oscillation wavelength proceeds by calculating the 
time average of probabilities over detection time, i.e., 
defining p a {L) ~ J Q T dt\A a (t, L)\ 2 , where T is the in- 
tegration time for the experiment. However, a summa- 
tion of probabilities over detection time is not justified 
by the rules of quantum theory. As the two-slit experi- 
ment indicates, alternatives that cannot be distinguished 
by measurement — in this case, detections of a flavor a 
at different moments of time — are to be summed over at 
the level of amplitudes. 

The summation over amplitudes, defined at different 
moments of time, is rather intricate because there is no 
time operator in quantum theory. Hence the methods 
that apply to other observables are not directly appli- 
cable here. The method we present in this letter has 
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the following advantages: (i) it fully implements the ba- 
sic principle, that alternatives that are not distinguished 
by measurement are to be summed at the level of am- 
plitudes; (ii) it incorporates a genuinely quantum treat- 
ment of detection time, according to the quantum theory 
of measurement [f|; (iii) it employs a general formula, 
Eq. ([3]), that is derived using only the rules of quan- 
tum theory. Eq. ([3]) is valid for any setup in which the 
detection time is not measurable, not only to particle 
oscillations. The end result is an expression for the de- 
tection probability as a function of distance; it leads to 
the non-standard formula for the oscillation wave-number 
kij = (rnf — rrij)/p that was mentioned previously. 

In what follows, we describe in more detail our ap- 
proach and the derivation of the results above. We em- 
ploy a method developed in Q, for the construction of 
probabilities from amplitudes that are defined at different 
moments of time. It has been applied to various prob- 
lems, such as probabilities for time-of-arrival, tunneling- 
time and for non-exponential decays. This method con- 
tains ideas from the decoherent histories approach to 
quantum mechanics Q and it has many similarities to 
the Srinivas-Davies photo-detection theory Q. 

The first step in our study is the derivation for a gen- 
eral formula for the probability of detection outcomes, 
when the time of detection is not observable. (The pre- 
sentation here is simplified: for elaboration on the finer 
points, see Ref. Q.) Let H be the Hilbert space of 
a quantum system, and H the Hamiltonian operator. In 
order to describe an event — such as a particle detection — 
we introduce a projection operator P that corresponds to 
states accessible only if the event has taken place. For 
example, if we detect a particle by monitoring a specific 
decay process, then P is a projector onto the states of all 
product particles. 

The Hilbert space H splits into two subspaces corre- 
sponding to the projectors P and Q = 1 — P. Let P\ be 
the projection operators corresponding to different values 
A of an observable that can be measured only if a detec- 
tion event has occurred. For example, P\ may correspond 
to a coarse-grained position variable for a product par- 
ticle. The set of projectors P\ is exclusive (P\P\' = 0, 
if A ^ A') and exhaustive, provided a detection has oc- 
curred; i.e., J2\P* = P- 

Starting from the initial state of a quantum system 
\ipo), w e construct the amplitude |?/>; t, A), at a final time 
T, that arises if a detection event took place within 
the interval [t, t + St) C [0, T] and the outcome of the 
measurement was A. By assumption, no detection took 
place before time t. Hence, we evolve Vo) with St — 
\iniN^, oc (Qe~ , ' Ht / N Q) N , i.e., with the restricted propa- 
gator in the eigenspace of the projector Q to no-detection 
states. In the interval [t, t + 5t] the state evolves accord- 
ing to the full unitary evolution e~ lHSt . Then, we act by 
a projection P\ that selects the measurement outcome 



A. Next, the state evolves unitarily until time T because 
there is no other event that needs to be taken into ac- 
count. At the limit of small St, the successive operations 
above yield \ip;t, A) = -iSt e~ lk ^ T '^P x HS t \^). 

We must emphasize here an important physical distinc- 
tion on the role of time pertaining to particle oscillations. 
The time of detection t is not identical to the evolution 
parameter of Schrodinger' s equation. Instead, it is a 
dynamical variable that determines the moment that a 
physical event has taken place • The construction of 
the amplitude \ip; t, A) above highlights this distinction: 
the detection time t is distinct from the time T at which 
the amplitude is evaluated. 

Furthermore, the amplitude \ip; t, A) is proportional to 
St, hence it defines a density with respect to time. The 
integration of \ip;t,X) over t is, therefore, well-defined in 
this scheme. In contrast, quantum theory provides no 
natural definition of integration of single-time probabil- 
ities over time Q- The amplitude \ip', X) corresponding 
to a detection with value A, at any (unspecified) time 
t e [0,T], is 

|^;A> = -i/ dte-^-^PxHSM-Cx^o). (2) 
Jo 

Hence, the probability p(X) that a detection with out- 
come A occurred at some time in [0, T] is 

= f dt f dt' Tr(e^ (t -^P x HSlp S t ,HP x ), (3) 
Jo Jo 

where po = \ipo)(4>o\- The probability measure Eq. Q is 
positive, linear, and normalized when the probability of 
no-detection 1 — L dtp(X) is included. 

If the Hamiltonian is of the form H = Hq + Hj, where 
[Ho, Q] = and Hj a small perturbation, then, to leading 
order in the perturbation, 

p(X) = f dt [ dt' TriUt^tPxHrUtpoUlHiPx), (4) 
Jo Jo 

where U t = e~ i6ot . 

Eqs. ((3HJ) are general operator expressions valid for 
any system. Next we apply them to the case of particle 
oscillations. Since the oscillating particles are detected 
by means of their decay products, the precise treatment 
of the detection process involves the use of quantum field 
theory. We will not specify the type of oscillating par- 
ticles so that our results will be valid for both neutrino 
and neutral boson oscillations. 

To this end, let us denote the oscillating particles as 
A and consider that their detection involves the process 
A + B m — » D n , where B m and D n are particles labeled 
by indices m and n. The Hilbert space H of the total 
system is a tensor product Ha ® Ti. r , where Ha is the 
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(bosonic or fermionic) Fock space J-(H\a)- The single- 
particle Hilbert space H\a is a direct sum (BiH-i of mass 
eigenspaces Hi- H r is the Hilbert space for the degrees 
of freedom corresponding to the B m and D n particles. 
It is decomposed as Ho © H pro d- Ho is the subspace of 
states prior to the decay A + B m — > D n and H pro( i is the 
subspace corresponding to states of the decay products 
[l2| . Since we assume that the measurements are carried 
out to the product particles, then the projection operators 
Pa in Eq. (J3j) are of the form 1 <S> tl\, where tlx projects 
into a subspace of H pro d- 

We assume a Hamiltonian of the form H = Ha ® 1 + 
l®H r + Hi, where Ha is the Hamiltonian for the A parti- 
cles, H r is the Hamiltonian for the B m and D n particles, 
and Hi is the interaction Hamiltonian. For simplicity, we 
assume that any particles B m that are present prior to 
detection are stationary. Hence, the restriction of H r to 
Ho is a constant, which can be taken equal to zero. The 

: equals eth+^2 n {\ 



We next consider the measurement of position X, of 
one of the product particles, with accuracy of order 5. 
The operators Ii\ in ((7]) can be substituted by a Gaussian 
approximate projector for position 



fix = / d 3 XV 



|X')(X'|®1, 



(8) 



restriction of H r on 7i 



prod 



MdJ, where e t h = J2n M D n ~J2 m M B m is the threshold 
energy of the A particles for the process A + B m — » D n . 
In the above, Ms m and Md„ are the masses of the parti- 
cles B m and D n respectively and p„ are the momentum 
operators for the D n particles. 

We consider an effective interaction Hamiltonian 



where the tensor product with unity refers to the remain- 
ing degrees of freedom in H pro d- 

The vector j+(x)|^ ) refers to the state of the prod- 
uct particles for decays that have taken place in a neigh- 
borhood of x. Hence, the operators Hx determine the 
locus of the decay event within an accuracy of order S, if 
ri(X)J+(x)|0 o ) ~ 0, for |x-X| » 5. Heuristically, the 
condition above would be satisfied for macroscopic values 
of 5 much larger than any length parameters character- 
izing the interaction. At this level of coarse-graining, the 
localization of one product particle at the decay time es- 
sentially determines the localization of all other product 
particles. In general, the localization scale S is macro- 
scopic, but it has to be much smaller than the scale 
of variation in the wave function of the A particles, or 
else no particle oscillations would be observable. Hence, 
to a first approximation, we can substitute ip a (x,t) by 
tp a (X.,t) in Eq. (J6j> . Taking these considerations into 
account, Eq. ((7|) becomes 



Hi = 



d 3 3 



b t (x)U at J+(x)+b\(x)U*J-(x) , (5) i?«(x,x',<-i')-^ 3 (X,x)<5 3 (X,x')e"-( t - t ')F(i'-i),(9) 



where bi , b\ are annihilation and creator operators on 
Ha, i labels mass eigenstates, (x) are current oper- 
ators of flavor a defined on H r , and Ui a is the mixing 
matrix (l3j|. The current operator involves products 
of annihilation operators for the B particles and creation 
operators for the D particles. Since no A particles are 
created during the detection process, the initial state \4>o) 
in H r must satisfy j~(x)|</>o) = 0. Note that if the de- 
tection involves a scattering process rather than a decay, 
an interaction Hamiltonian quadratic to the field of the 
A particles should be used, instead of (|5|). 

Let the initial state on Ha be a single-particle state 
IV'o) =J2iI d 3 xbl(x)tjj i0 (x.)\0) a, where \0) A is the vac- 
uum of Ha- Since [H ,Q] = 0, Eq. Q applies. Hence, 
the probability that a decay through flavor a has hap- 
pened at some time in [0, T] and that the value A for an 
observable of the product particles has been found equals 

Pa( x ) = J2f o dt f Q dt ' J d 3 xd 3 x'i(}](x',t')i;i(x,t) 

xU ai U* a R%(x,x',t-t'),(6) 

where ipi (x, t) is the evolution of V^o ( x ) under the Hamil- 
tonian for a single A-particle and 



where K > is a constant. The function F(s) is 
obtained from the propagator of H r and it equals 

[I J d a p c v n ; n runs over all 

product particles D n . If M Dn S » 1, the saddle-point 
approximation applies and 



F(s) = n 



Mr 



2m{s-iM D J 2 /2) 



3/2 



(10) 



Eqs. ([6]) and ([9]) define a probability density for the 
particle's position x at scales much larger than 8 

p a {x) = K / dt / dt'A a {t l x)A* ol (t' > x)e- i ^ t '^F(t'-t)(ll) 



where A a (t,x) — ^2nU a itpi{t, x) is the probability am- 
plitude corresponding to the flavor a. The probability 
in Eq. (|11[) can be explicitly computed for any type of 
oscillating particle, provided we specify its initial state. 

To this end, we consider the evolution of an initial 
wave-packet of A particles. We assume that the solid an- 
gle connecting the production region and the detection 
region is very small, so that only particles with momen- 
tum along the axis that connects the two regions are de- 
tected. Let the initial state of the A particles be a flavor 
superposition of different Gaussian mass-eigenstates 



x',t- 



-t') = (0o|i-(x')n Ae -^ (t '" t) n A j+(x)|0 o ).(7) 



tp i0 {x) = U* 0i 



1 



(ira 2 ) 1 / 4 



T-\-iptx 



(12) 
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where a is the spread of the initial wave-packet and (3 
the initial flavor. Then, 



/4 



A a (t,x) =J2u* 0i U ai (47ra 2 ) 1 



dp ( P - Pi ) 2 +ipx-iE t (p)t-r< ( P )t 
2tt 



(13) 



where Ej (p) = */ mf + p 2 and 1^ (p) are the decay rates in 
the different mass eigenspaces. We employ a commonly 
used approximation: we expand Ei(p) to first order in 
p — pi and Ti(p) to zero-th order in p — pi, i.e., 



Ei(p) = Ei+ Vi(p-pi); 



T i (p)=T l 



(14) 



where Ei = Ei(j>i), Vi = Pi/Ei, and I\ = Ti(pi) \U 
Then, 



- (x "'^ — MpiX-iEit-Fit 



(trt 2 ) 1 ^ 



(15) 



Substituting (fl~5f in Eq. (fTTjl and letting T — > oo, we 
obtain the probability density for the detection of flavor 
a at x = L 

Pa (L) = K ' (j2 S ^ L + E 2e " ( " 

xi? e (T„- e lfc - L )) , (16) 



where if ' > is a redefined normalization constant, Si ~ 
[/*-l7 m | 2 , T 4i ~ U^UaiUpjU*] (their precise form is not 
necessary for the arguments) and 



fey — 



Ei — e th Ej — e th 



-(Pi-Pj). (17) 



If eth/Ei << 1 then k i3 = m 2 /p l - m 2 /pj, i.e., for pi = 
p+0(mi), kij is twice the standard result (m 2 — m 2 )/(2p). 

It is important to emphasize that this result does not 
depend on any ad-hoc assumptions about the initial state. 
It follows from Eq. ([3]), which is valid for a large class 
of measurements, and it provides the probabilities for 
alternatives when the detection time is not measured. 
Moreover, Eq. (jTTJ) is insensitive to the approximations 
we made in this paper. It holds for any calculation, in 
which the summation over detection time is performed 
at the level of amplitudes. This is the reason why it it is 
discerned even in the simplified plane wave description. 

If we do not assume that eth/Ei « 1, Eq. (fTT|) be- 
comes = (l — ^)(m 2 —m 2 )/E, in the ultra-relativistic 
limit. The oscillation wavelength carries a strong depen- 
dence on the threshold energy. This dependence does 
not arise at all in the standard treatment and it is not 
an artifact of any approximation. In particular, the pres- 
ence of an energy threshold is equivalently described as 



a constant potential V — —eth, in which the particle 
propagates prior to detection. The incorporation of this 
potential to the Hamiltonian would lead to Eq. (fT6|) even 
in the simplified plane wave description. 

As a final remark, we note that the kernel F[t — t') in 
Eq. (jTTJ) turns out to cause significant suppression of in- 
terferences in detection time for the amplitude (I15p . One 
might, however, inquire whether it is possible that in a 
different model one might obtain F(t — t') ~ S(t — t'), 
so that p a {L) ~ J Q T dt\A(t, L)\ 2 \ the standard expression 
for the oscillation wavelength would then follow instead 
of (|17[) . However, this would require the presence of an 
unusual and highly efficient decohering mechanism for 
the detection time that would be effective in timescales 
shorter than the shortest one appearing in (|15j) . namely, 
EJ . Moreover, such a mechanism would have to be ex- 
trinsic to the physics of particle oscillations, and it would 
have to be postulated ad hoc. 
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